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Abstract 

We present a calculation of the leading order QCD fragmentation func¬ 
tions for gluons to split into spin-triplet D-wave quarkonia. We apply 
them to evaluate the gluon fragmentation contributions to inclusive 3 Dj 
quarkonium production at large transverse momentum processes like the 
Tevatron and find that the D-wave quarkonia, especially the charmo- 
nium 2 state, could be observed through color-octet mechanism with 
present luminosity. Since there are distinctively large gaps between the 
contributions of two different (he, color-singlet and color-octet) quarko¬ 
nium production mechanisms, our results may stand as a unique test to 
NRQCD color-octet quarkonium production mechanism. 
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I. INTRODUCTION 


The production and decays of quarkonium bound states have been under active studies in 
experiment ever since the first charmonium 1 state, the J/ip, was found twenty years before. 
The study of properties of the bound states of heavy quarks has provided a wealth of information 
on both the properties of heavy quarks (c, b ) themselves and quantum chromodynamics since 
it stands on the very border between perturbative and nonperturbative domains. Recently, the 
observations of the Collider Detector at Fcrmilab (CDF) on the prompt charmonium production 
has greatly stimulated progress in theoretical studies of quarkonium physics. 

The conventional wisdom was that the dominant contributions to quarkonium production 
cross section at large transverse momentum (Pt) hi pp collisions come from the QCD leading 
order diagrams, i.e. the so called parton fusion processes. However, these calculations for 
prompt J/ip (ip') production did not reproduce all the aspects of the available data [||. It 
was pointed out by Braaten and Yuan [|3] in 1993 that the dominant production mechanism at 
sufficiently large Pt is the fragmentation of a parton produced with large transverse momentum, 
while formally this is of higher order in the strong coupling constant a s . Unfortunately, even 
after including the fragmentation contributions, the predictions for the ip' production rate still 
falloff far below the data [(4|. This large discrepency between theory and experiment has called in 
question the simple color-singlet model description for quarkonium || and suggests that a new 
paradigm for treating heavy quark-antiquark bound systems that go beyond the color-singlet 
model might play an important role in the production of quarkonium at large Pt- 

To this end, a factorization formalism has recently been performed by Bodwin, Braaten, 
and Lepage [|6|] in the context of nonrelativistic quantum chromo dynamics (NRQCD), which 
provides a new framework to calculate the inclusive production and decay rates of quarkonia. 
In this approach, the calculations are organizeded in powers of v, the average velocity of the 
heavy quark (antiquark) in the meson rest frame, and in a s , the strong coupling constant. 
In NRQCD, quarkonium is not solely regarded as simply a quark-antiquark pair but rather a 
superposition of Fock states. The general Fock state expansion starts as 

\H(nJ PC ) > = 0(1)\QQ( 2S+1 Lj,1) > (1) 

+ 0(v)\QQ( 2S+1 (L±l)j,,8)g> 

+ 0(v 2 )\QQC s+1 Lj,8)gg>+--- 

H- 
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where the angular momentum quantum numbers of the QQ pairs within various Fock com¬ 
ponents are indicated in spectroscopic notation inside the brackets with color configuration of 
either 1 or 8. 

The breakdown of color singlet model stems from its overlook of high Fock components 
contributions to quarkonium production cross sections. The color-octet term in the gluon frag¬ 
mentation to J/ip(ip') has been considered by Braaten and Fleming Q to explain the 
surplus problems discovered by CDF. Taking < (D$^( 3 Si) > and < Og ( 3 Si) > as input pa¬ 
rameters, the CDF surplus problems for J/0 and 0' can be explained as the contributions of 
color-octet terms due to gluon fragmentation. 

Even though the color-octet mechanism has gained some successes in describing the pro¬ 
duction and decays of heavy quark bound systems ili it still needs more effort to go 
before finally setting its position and role in heavy quarkonium physics. Therefore the most 
urgent task among others needs to do now is to confirm and identify the color-octet quarkonium 
signals. 

While the Erst charmonium state, the «//-0, has been found over twenty years, D -wave states, 
given the limited experimental data, have received less attention. However, this situation may 
be changed in both experimental and theoretical investigations. Experimentally, there are 
hopes of observing charmonium D -wave states in addition to the -0(3770), in a high-statistic 
exclusive charmonium production experiment [9|, and bb D -wave states in T radiative decays 

0 . 


Recently, there is some clue for the D -wave 2 charmonium state in £705 300 GeV 7r - and 


proton- Li interaction experiment |lT|. In this experiment there is an abnormal phenomenon 
that in the J/07r + 7r _ mass spectrum, two peaks at 0(3686) and at 3.836 GeV (given to be 
the 2 state) are observed and they have almost the same height. Obviously, this situation 
is difficult to explain based upon the color-singlet model. However, it might be explained 
with the NRQCD analysis. Of course, at energies in fixed target experiments like £705, the 
color-octet gluon fragmentation dominance may or may not be the case. Moreover, the strong 
signal of J/07r + 7r _ at 3.836GeE observed by £705 is now questioned by other experiments 
0 - Nevertheless, if the £705 result is confirmed (even with a smaller rate, say, by a factor 
of 3, for the signal at 3.836GeE), the color-octet gluon fragmentation will perhaps provide a 
quite unique explanation for the D -wave charmonium production. It does remind us that in the 
NRQCD approach, as discussed in Ref. [I3|, the production rates of /7-wave heavy quarkonium 
states may be as large as that of F-wave states as long as the color-octet gluon production 
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mechanism dominantes. 

Study shows that both LEP and Tevatron, especially the latter, are suitable grounds to 
find the D-wave quarkonia [|T3| and to test the color-octet signals |T4|]. In this paper we find 
that the divergences of the contributions between color-singlet and color-octet mechanisms in 
quarkonium 3 Dj production is enormous. The rest of the paper is arranged as follows: In Sec. 
II, we describe the formalism, give out the fragmentation functions of g —> 3 Dj to leading order 
in a s and further calculate the fragmentation probabilities of the gluon to D-wave quarkonium 
states. In Sec. Ill, we apply the fragmentation functions to evaluate the D-wave quarkonium 
production rates at Tevatron and close with some thoughts and discussions. 


II. FORMALISM 


Fragmentation is the formation of a hadron within a jet produced by a parton (quark, anti¬ 
quark or gluon) with large transverse momentum. It is a useful concept because the probability 
for the formation of hadron within a jet is independent of the process that produces the parton 
that initiates the jet. By now, the fragmentation functions of quark and gluon splitting to S- 
and P-wave heavy quark bound states have been calculated [[Uj |H| |T6| ]I7[. The calculations 


of fragmentation functions of quark to D-wave states |18| and gluon to spin-singlet D-wave 
state 1 D 2 have also been accomplished 0. However, the study of gluon fragmentation to 
color-singlet 3 Dj states, for its complexity, is still left behind. Here we realize this goal. 

In hard process, as Fig. 1(a), the most important kinematic region for a virtual gluon split 
with large Pt is that the gluon is nearly on its massshell. Therefore, we can estimate the decay 
widths and the branching ratios by the following way |20[ . 

The decay widths of a virtual quark Q* to color-singlet quarkonium state 3 Dj by gluon 
fragmentation can be evaluated via 


S 

r(Q* - Qg*;g* - 3 Dj gg ) = J d g 2 r(Q* - Q g *(»)) ■ P ( g * -■ 3 d, ; gg ), 


( 2 ) 


where s is the invariant mass squared of Q*\ p is the virtuality of the gluon, and its minimum 
value squared gP n = 12 mg corresponding to the infrared cutoff as discussed below; P is the 
decay distribution defined as 


P(g*^AX) = —T{g*^AX). 

7T/P 


( 3 ) 
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The general covariant procedure for calculating the production and decay rates of heavy 
quark bound states may start from the Bethe-Salpeter(BS) amplitudes in the nonrelativistic 
limit. At leading order in a s , the amplitudes for g* ~^ 3 Dj gg processes are 

A = f -0^ Tr {O( p ,q)x( p ,q)}- (4) 

Here y(P, q ) is the BS wave function of the bound states with relative momentum q between 
the heavy quarks, while 0(P , q) represents the rest of the matrix elements depicted in Fig.l(a), 

- fti+ ftz+ ft ~ 2 m Q - lfti+ ft 2 - ft - 2 m Q 

° {<! > = ~4 { & -(h -k)-k 2 <■ * } 

+ five permutations of Ay, k 2 , —k and e±, e 2 , e . (5) 


Here the Aq ; k 2l k and ei, e 2 , e stand for the momenta and polarization vectors of the two 
outgoing final gluons and the splitting gluon. Coupling constants and color matrices have been 
supressed and contribute a factor 


E(4 5 9, s rUr“r'T'}) 2 = A,, 6 

a.b.c V 3 lo 


( 6 ) 


to the production rates. 

Under the instantaneous approximation with the negative energy projectors being neglected, 
the BS wave function y(P, q) may be expressed as 


x(p, q) = 


Po — Ei ~ E 2 


MP,q)- 


(7) 


2vr (p w - E 1 )(p 20 - E 2 ) 

Here Po is th e time component of the four momentum of the bound state; p 10 and p - 2 o are the 
time components of the momenta of quark and antiquark inside the meson, and Pi, E 2 are their 
kinetic energies. From the standard BS wave functions in the approximation that the negative 
energy projectors are omitted, the vector meson wave function can be projected out as : 


1 


$(P,q) = jj ^( JM |l^Pm)A+(pl)7o ft(M+ f) 7o A 2 _(p 2 )^ m (P,g), 


( 8 ) 


S z m 


where e is the polarization vector associated with the spin-triplet states. A^(pj) and A^(jq>) 
are positive energy projection operators of quark and antiquark . 


Al /^s Ei + 7o7 • Pi + ^i7o a2 /^v E 2 - 7 q 7 • p 2 - m 2 7o 
A+(pi) =-TP-1 A _(p 2 ) = 


(9) 


2Pi ’ 7 2P 2 

After taking the nonrelativistic approximation the bound state wave function may be further 
reduced. For D-wave quarkonium production and decay, the first nonzero term is proportional 
to the second order of the amplitude expansion in powers of q/M: 
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( 10 ) 


A(P, q) = A(P, 0) + q c 


9A(P , q ) 
dq a 


=o + 


19=0 


1 d 2 A{P , q) 

2 q a q q Q qa Q q p 


q =0 + ' ' ' • 


After integrating q a qg over d 4 g, the 3 Dj polarization tensor is related to its nonrelativistic 
wavefnnction by 


/ 


d 3 q 

(2tt) 3 


q a qpil’ 2 m(P, q) 



( 11 ) 


where the polarization tensor’s label m ranges over the helicity levels of the L — 2 meson. 
For spin-singlet case is identified with e^g , while for the spin-triplet case, using explicit 
Clebsch-Gorden coefficients, we have the following spin-orbit momentum coupling forms |j21|j, 


Y,( 1J z\ 1S z 2m )e < $e { * ) 

S z m 


Y J {^Jz\lS z 2m)e { ^e^ ) 

S z m 

( 3 . 7 Z | 13*2 m)e { ^e^ 


S z m 


\ ^ ] 1/2 r/ PaPp\(J z ) , / PPPp\AJz) 

20 ^ +( 9* W )£ « 


Q 

_ PaP d \ e (Jz) 1 


Q 


2^/6m Q 


l e acr CrPpcr'P 9 1 £/3cr e Taper' V 9 ) j 


,(J») 

t a/3p- 


( 12 ) 

(13) 

(14) 


Using Eqs. (|T2D-(|T4|) listed above, the amplitudes of Eq. ([4]) may be simplified and the averaged 
squared amplitudes may be obtained when suming up all polarizations of both the meson and 
gluons. Because the results are lengthy, it is too tedious to write them all here. For the 
convenience of reference, we just give the expression for the 3 D\ state in the Appendix. Then, 
we have 


r(s - —* 3 Dj gg) = JdjgdjgX\A 2 . 


(15) 


where the kinematic variables are defined X\ = and = ^Ap-. Furthermore, from the 
Eq.(jKJ) we can get the expressions of decay distributions P(g* —> 3 Dj gg). With them the 
fragmentation functions can be calculated straightforward 

dT(Q* —> 3 Dj gg Q)/dz 


D g *^3 Dj (z,2m Q ,s) = 


(16) 


r(Q* - Qg) 

where z = = 2 — x± — rr 2 - At high energy limit, the interaction energy s goes up to infinity, 

then the definition of z here is identical with that in Ref. |]T7[ multiplied by a factor of two and 
the fragmentation functions decouple from any specific gluon splitting processes, which just 
reflects the universal spirit of fragmentation. The fragmetation function of Eq. (JlGi) is evaluated 
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at the renormalization scale 2 rriQ, which corresponds to the minimum value of the invariant 
mass of the virtual gluon. In Fig.2 we display the variation curves of D g *^3 Dj (z,2m c ) versus 
z. After integrating over variable z, the fragmentation probabilities then read as 

r(Q* Dj gg Q ) 


Pn 


■» 3 D, ~ 


(17) 


r (Q* - Qg ) 

Studies show [17] that the above method in extracting the gluon fragmentaiton probabilites are 
equivalent to the method developed in Ref. 

The calculation of color-octet fragmentation functions in g* —> 3 Dj( 3 S 1; 8) processes, as 
shown in Fig.l(b), is trivial. They may be obtained directly from color-octet g* —> J/ip^S i,8) 
process 


D g *^3 Dj (z,2m Q ) = 


na s (2m Q ) 
24 m 2 Q 


5(1-z)<(D 8 d J( 3 S 1 )> . 


Therefore, the fragmentation probabilities are expressed as: 


Pn 


* 3 Dj ~ 


^#^4 < ol D >es o >. 


24 m 2 Q 


(18) 


(19) 


III. RESULTS AND DISCUSSIONS 


From Eq.([17]) and (|l^) we can estimate the quarkonium 3 Dj production rates at the Teva- 
tron. The color-singlet sector may be factorized into long distance and short distance terms. 
The former is, to leading order in v 2 , proportional to the second derivative of the radial wave 
function at the origin, which may be determined from potential model calculations |^2|. The 
latter can be calculated from perturbative QCD, and it involves the infrared divergence asso¬ 
ciated with a soft gluon in the final states. In the numerical computation, we impose a lower 
cutoff A on the energies of either gluons in the quarkonium rest frame. As discussed in Ref. 

T7j| , we choose A = tuq to avoid large logarithms and the cutoff dependence of the color-singlet 
terms is cancelled by the A dependence of the nonperturbative matrix elements < C ) 8 i:,j ( 3 S'i) > 
of the corresponding color-octet terms. 

The gluon fragmentation contributions to the production of quarkonium 3 Dj states at large 
transverse momentum in any high energy process can be approximately obtained by multi¬ 
plying the cross section for producing gluons with transverse momentum larger than 2 m c by 
appropriate fragmentaion probabilities 0. Using |T5[ [|22 


m c = 1.5 GeV,mb = 4.9 GeV, a s (2m c ) = 0.26, a s (2m b ) = 0.19, 
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( 20 ) 


We obtain 


R" cg) (0)\ 2 = 0.015 GeV 7 , | (0) | 2 = 0.637 GeV 7 . 


]jG) _ x a ^ in-8 nh) 


+ 3D l( cc) = 5-6 x 10- 8 , D^ 3D2{c5) 


D 


(i) 

g*^ 3 D 3 (cc) 


= 2.2 X 10 


-7 


D 


(i) 

g*—> 3 D\ ( bb ) 


Dg}_>3D 2 (bb) 1-4 x 10 , -Dg*^3 D3 ( 6fe ) 


: 3.1 x 10“ 7 , 
2.5 x 10" 10 , 
9.9 x 10" 10 . 


( 21 ) 


For gluon fragmentation color-octet processes, the fragmentation probabilities are propor¬ 
tional to the nonperturbative matrix elements < e>8 Dj ( 3 Si) > which have not been extracted 
out from experimental data, nor from the Lattice QCD calculations. Based upon the NRQCD 
velocity scaling rules and the experimental clues discussed above, here we tentatively assume 

in 

< e>g D2(ca) ( 3 Si) >«< of ( 3 Si) >= 4.6 x 10” 3 GeV 3 (22) 


(see Ref. 0) and further extend this relation to the bb system 0 

< Os D2(b " b) ( 3 Si) >~< 0j'( 3 Si) >= 4.1 x 10 -3 GeV 3 . (23) 


The supposed relations (|22|) and (^3j) certainly possess uncertainties to some extent, however 
from the calculated results below we are confident that it will not destroy the major conclusion 
of this paper. From the approximate heavy quark spin symetry relation, we have 

< Os Dl ( 3 Si) >~ l < Ol D2 ( 3 S 1 ) >^ < Ol^CSi) > (24) 

5 7 

for both bb and cc systems. 

Using Eqs.(|l9D, fl20|), and (j22l)-(j24p, we readily have 


£><? 

£><? 

9 ~ 


3 D\(cc) 


3 D 3 (cc) 



3 D 2 (bb) 


4.2 X 10“ 5 , D {s) 


’ g*^ 3 D 2 (cc) 


9.7 x 1(T 5 , D [ °J 3n 

7 g —>*D i 


(8) 


4.2 x hr 6 , D ( l 

* n* 


(bb) 

(8) 

g*^3 D3 ( b b) 


7.0 x 10~ 5 , 
2.5 x 10~ 6 , 
5.9 x hr 6 . 


(25) 


Comparing the above results (|25|) with ([H|), we come to an anticipated conclusion that at 
the Tevatron the gluon fragmentation probabilities through color-octet intermediates to spin- 
triplet D-wave charmonium and bottomnium states are over 2 ~ 4 orders of magnitude larger 
than that of color-singlet processes. As a result, the production rates of 3 Dj states are about 
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the same amount as 0' and Y(2 S) production rates. Compared with the 0' production at 
the Tevatron, the gluon fragmentation color-octet process plays an even more important role 
in the 3 Dj quarkonium production, and it also gives production probabilities larger than the 
quark fragmentation process 


Among the three triplet states of D -wave charmonium, 3 D- 2 is the most promising candidate 
to discover firstly. Its mass falls in the range of 3.810 ~ 3.840 GeV in the potential model 
calculation m, that is above the DD threshold but below the DD* threshold. However the 
parity conservation forbids it decaying into DD. It, therefore, is a narrow resonance. Its main 
decay modes are expected to be, 


3 D 2 -> J/07T7T, 3 D 2 


► 3 P j7 (J = 1,2), 3 D 2 


3 g. 


(26) 


We can estimate the hadronic transition rate of 3 D 2 —> J/07r + 7r _ from the Mark III data for 
0(3770) —> J/07T + 7r _ |24] and the QCD multipole expansion theory |25| p6| . The Mark III 
data give [p| T(0(377O) —> J/07r + 7r _ ) = (37 ± 17 ± 8) keV or (55 ± 23 ± 11) keV (see 
also Ref. pi). Because the S — D mixing angle for 0(3770) and 0(3686) is expected to be 


small (say, —10°, see Ref. f27j for the reasoning), the observed 0(3770) —> J/07r + 7r“ transition 
should dominantly come from the 3 D\ —> transition, which is also compatible with 

the multipole expansion estimate p||. Then using the relation [p| 


dT( 3 D 2 = dT( 3 D 1 —0 S\2tv) 

and taking the average value of the T(0(377O) —> J/07r + 7r _ ) from the Mark III data, we may 
have 


T( 3 D 2 -»■ J/07T+7T-) = r( 3 D! -»■ j/07T+7T-) « 46 keV. 


(27) 


For the El transition 3 D 2 —> 3 Pjj^J = 1,2), using the potential model with relativistic effects 


being considered [28], we find 


T( 3 A> -> Xdl) = 250 keV, T ( 3 D 2 ->• Xc 2 l) = 60 keV, 


(28) 


where the mass of 3 D 2 is set to be 3.84GeE. As for the 3 D 2 —> 3 g annihilation decay, an 


estimate gives 29 


T( 3 D 2 -> 3 g) = 12 keV 


(29) 


From (|27|) , (|28|) , and (|29|) , we find 
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( 30 ) 


r«( 3 B 2 )» r( 3 £> 2 j/i>i n) + r( 3 £> 2 — Xc a) + r ( 3 b 2 -> x*n) + r( 3 c 2 39) 

rs 390 keV, 


and 


B( 3 D 2 —■> J/'07T + 7T ) ~ 0.12. (31) 

Considering all the uncertainties this estimate is expected to hold within 50%. Compared (PH) 
with B{ip' —■> J/07r + 7r _ ) = 0.324 ± 0.026, the branching ratio of 3 D - 2 —> J/07r + 7r^ is only 
smaller by a factor of 3, and therefore the decay mode of 3 D 2 —> J/07r + 7r _ could be observable 

at Tevatron. 

The 3 Di cc state -0(3770) could also detected via 3 Di — > DD. The other states, including 
the 3 D 3 (cc ) and 3 Dj(bb) are perhaps difficult to detect for reasons of either more decay modes 
or smaller production rates. 

In conclusion, we have calculated the fragmentation functions and fragmentation proba¬ 
bilities of the gluon to 3 Dj charmonium and bottomonium states in both color-singlet and 
color-octet processes with certain numerical assumptions (e.g. Eq.(^)). The results can also 
be used in other hard gluon fragmentation processes because of the universality of the fragmen¬ 
tation functions. The study shows that, because charmonium 3 D 2 state may have a production 
rate as large as that of 0' at the Tevatron through color-octet production mechanism, the 
charmonium 3 D 2 state as a most promising candidate to discover should be observable at the 
Tevatron with present luminosity, even the assumption of Eq.(^) with an error of 10 times 
off the exact case. On the other hand, since the calculated results show that the color-singlet 
and the color-octet contributions diverge enormously, this will also present a crucial test for 
the color-octet mechanism, the 3 Dj bottomonium states may have less strong signals to be 
detected comparing with the 3 Dj charmonium states because of their small production rates. 
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APPENDIX 


The expression of the averaged squared amplitude of the process g* —> A Digg 


J2\A\ 2 = 25 a 2 ^\ R D (°)\ 2 


128 d £ fid { 

i =0 


2 15 7Tm 7 c 15(1 — d — xi) 5 (l — d — x 2 ) 5 ( a; i + x 2 ) 6 ’ 

where d = M 2 / g 2 and the functions f t are defined as 


( 32 ) 


/o = (a+ — l) 4 (x 2 — 1) 4 [Ti(x 2 + 17) + x^(4x 2 + 53x 2 — 177) + xf(6xl + 74x 2 — 420x 2 
+576) + x\{Ax\ + 74x 3 - 486x^ + 1216x 2 - 928) + xi(x^ + 53x^ - 420^ + 1216x^ 
-1600a; 2 + 768) + 17x\ - 177x2 + 576x 3 ~ 928x 2 + 768x 2 ~ 256], 
fi = {x\ — 1) 3 (x 2 — l) 3 [10xf + 4 x]’(15x 2 + 1) + x®(170x 2 + 39x 2 — 141) + xf (300x 3 
—133x2 - 6i6x 2 - 271) + x^(360x^ - 566x 3 - 803x 2 - 587x 2 + 3162) + x 3 (300x® 
—566x2 - 656x2 + 106x^ + 8336x 2 - 8696) + x?(170x^ - 133x| - 803x^ + 106x^ 
+9772x2 - 21480x 2 + 12720) + xi(60x^ + 39x® - 616x 3 - 587x* + 8336x 3 
—21480x2 + 24160x 2 - 9856) + 10x^ + 4x^ - 141x£ - 271x^ + 3162x^ - 8696x^ 
+12720x^- 9856x2 + 3136], 

f 2 = (xi - 1) 2 (x 2 - l) 2 [-50x? + 2x 3 (—125 x 2 + 14) + 2x[(-320x^ + 479x 2 - 97) + x?(1951 
+3593x2 - 5120x2 - 1120x 3 1951) + 2x^(-730x^ + 2819x 3 - 8038x 2 + 7736x 2 - 160) 
+x^(—1460x2 + 5950x2 - 24146x^ + 44649x^ - 11984x 2 - 17579) + 2x?(-560x^ 
+2819X® - 12073x2 + 31416x 3 - 22632x 2 - 24258x 2 + 26192) + x 2 (3593x® - 640x^ 
—16076x2 + 44649x2 - 45264x^ - 57330x^ + 146496x 2 - 75496) + 2xi(479x^ - 125x^ 
—2560x2 + 7736x2 - 5992x^ - 24258x^ + 73248x5; - 76552x 2 + 27952) - 50x|] + 28x| 
—194x2 + 1951X® - 320x2 - 17579x]] + 52384x5] - 75496x5] + 55904x 2 - 16832], 

/ 3 = 2(xi - 1)(x 2 - l)[67xf + x?(277x 2 - 251) + 3x?(328x5] - 618x 2 + 405) + xj(2780x5] 

—5699x 2 + 3156x2 - 513) + x?(5153x(] - 12461x5] - 1731x5] + 18247x 2 - 8577) 

+x 3 (6294x2 - 19575x2 ~ H940x| + 89313x5] - 82353x 2 + 15559) + x?(5153x£ 
—19575x2 - 16760x2 + 178137x5] - 272631x5] + 116247x 2 + 11894) + x?(2780x]] 
—12461X® - 11940x2 + 178137x2 - 401134x5] + 303362x5] + 20372x 2 - 79388) 
+x 2 (984x® - 5699x2 - 1731x5] + 89313x5] - 272631x5] + 303362x5] + 10876x5] 
-249684x2 + 125056) + xi (277x5] - 1854x5] + 3156x5] + 18247x5] - 82353x5] 

+116247x2 + 20372x2 - 249684x5] + 266656x 2 - 91040) + 67x]° - 251x5] + 1215x5] 
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-5134 - 85774 + 155594 + 118944 - 793884 + 1250564 - 91040a; 2 + 26080, 
u = 98a;} 1 + 2a;}°(629a; 2 - 640) + 4(68044 - 13203a; 2 + 6503) + 4(236124 

-635064 + 56515a; 2 - 17151) + 2x^(266374 - 901834 + 976724 - 32167a; 2 
-1868) + 24(396214 - 1671394 + 2059924 - 28914 - 135604a; 2 + 59500) 
+24(396214 - 2053994 + 2962354 + 1423714 - 7021844 + 558303a; 2 - 127486) 
+24(266374 - 1671394 + 2962354 + 2464774 - 14844084 + 18355804 
—839002a; 2 + 84843) + 24(118064 - 901834 + 2059924 + 1423714 - 14844084 
+26894584 - 19788244 + 403532a; 2 + 100128) + 24(34024 - 317534 + 976724 
—28914 - 7021844 + 18355804 - 19788244 + 6405464 + 380960a; 2 - 242436) 

+xi (12584° - 132034 + 565154 - 643344 - 2712084 + 11166064 
-16780044 + 8070644 + 7619204 - 1085712a; 2 + 369120) + 984 1 - 12804° 
+65034 - 171514 - 37364 + 1190004 - 2549724 + 1696864 + 2002564 
—4848724 + 369120a; 2 - 102720, 

/ 5 = 490a;} 0 + x?(4870x 2 - 4351) + 4(250964 - 42443a; 2 + 17526) + 24(389604 
-937614 + 63537a; 2 - 8870) + 24(756954 - 2409894 + 2076274 + 12406a; 2 
-54294) + 24(941864 - 3831334 + 4071234 + 1938184 - 526995a; 2 + 209916) 
+24(756954 - 3831334 + 5070364 + 4805024 - 17222864 + 1369912a; 2 
-320530) + 24(389604 - 240989a; 0 + 4071234 + 4805024 - 25150424 
+31943644 - 1563232a; 2 + 196056) + 24(125484 - 937614 + 2076274 
+1938184 - 17222864 + 3194364a; 0 - 25118044 + 648200a; 2 + 70608) + a;i(4870xjj 
-424434 + 1270744 + 24812a; 0 - 10539904 + 27398244 - 3126464 x\ 

+12964004 + 360384a; 2 - 330624) + 490x}° - 4351a; 0 + 175264 - 177404 
-1085884 + 4198324 - 6410604 + 3921124 + 1412164 ~ 330624a; 2 + 131712, 
f 6 = 9404 + 4( 9380a; 2 - 6619) + 24(237544 - 33091a; 2 + 8148) + 24(670624 

—1423034 + 55480a; 2 + 22478) + 24(1110644 - 3274854 + 1786204 + 209904a; 2 
-168145) + 24(1110644 - 4296394 + 3100244 + 6671504 - 1059801a; 2 + 386365) 
+24(670624 - 3274854 + 3100244 + 9684084 - 24329024 + 1859180a; 2 
-431760) + 24(237544 - 142303a; 0 + 1786204 + 6671504 - 24329024 
+29810544 - 1469776a; 2 + 192024) + 2a;i(46904 - 330914 + 554804 + 2099044 
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-10598014 + 18591804 - 14697764 + 388080x2 + 45024) + 9404 - 66194 
+162964 + 449564 - 3362904 + 7727304 - 8635204 + 3840484 
+90048x 2 - 104832, 

f 7 = 2[410x® + 3x^(1849x2 - 804) + 2x) i ( 14361x5] - 15578x 2 - 1348) + x?(71493x5], 
—131020x2 - 7327x2 + 72719) + x^(95816x5) - 253812x5] - 3044x5] + 426335x 2 
-254736) + x 5] (71493x5) - 253812x5] + 582x5] + 950322x5] - 1194424x 2 + 406920) 
+2x^(14361x2 - 65510x2 - 1522x5] + 475161x5] - 951640x5] + 691660x 2 - 157488) 

+xi (5547x2 - 31156x2 - 7327x5] + 426335x5] - 1194424x5] + 1383320x5] 

-661824x2 + 77952) + 410x5] - 2412x5] - 2696x5) + 72719x5] - 254736x5] 

+406920x2 - 314976x2 + 77952x 2 + 19392], 
fs = 454x1 + 2x^(5231x2 - 879) + x^(49146x5] - 37439x 2 - 29865) + x^(98210x5] 

—138246x2 - 144901x2 + 193291) + 2x^(49105x5] - 102277x5] - 148785x5] + 435590x 2 
-221376) + 2x^(24573x2 - 69123x5] - 148785x5] + 686625x5] - 744064x 2 + 240264) 
+Xi( 10462x2 - 37439x2 - 144901x5] + 871180x5] - 1488128x5] + 1017888x 2 - 221760) 
+454x2 - 1758X® - 29865x5] + 193291x5] - 442752x5] + 480528x5; - 221760x 2 + 13440, 
fg = 550x® + x®(8230x 2 + 637) + x^(28258x5] - 6151x 2 - 38526) + 2x^(20578x5] 

—10371x2 - 81952x2 + 71372) + 2x^(14129x5] - 10371x5] - 126946x5] + 234980x 2 
-106504) + xi (8230x5) - 6151x5] - 163904x5] + 469960x5] - 451168x 2 + 140480) 
+550x5] + 637x2 - 38526x2 + 142744x5] - 213008x5] + 140480x 2 - 25280, 

/io = 572x^ + x^(4036x 2 + 2537) + 4x^(2288x5] + 2455x 2 - 6348) + 2x^(4576x5] 

+7379x5] - 40504x2 + 27420) + 4xi( 1009x5] + 2455x5] - 20252x5] + 28844x 2 - 12392) 
+572x5] + 2537x2 ~ 25392x5] + 54840x5] - 49568x 2 + 14144, 
fn = 4[51x? + 2x?(101x 2 + 245) + 2x^151x5] + 751x 2 - 928) + 2 xi(101x5] 

+751x^ - 1920x2 + 1152) + 51x5] + 490x5] - 1856x5] + 2304x 2 - 976], 
f V2 = 8 [51x2 + 6 xi(17x 2 - 14) + 51x5] - 84x 2 + 56], 
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Figure Captions 


Fig.l. Virtual gluon fragmentation processes (a) gluon fragments to 3 Dj via color-singlet pro¬ 
cess, (b) gluon fragments to 3 Dj via color-octet process. 

Fig.2 The variation of charmonium fragmetation functions Dr g _>3 Dj \(z,2m c ) versus z. 
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